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PACS 36 . 20 . -r - Macromolecules and polymers 
PACS 64 . 75 . Yz - Self-assembly 

PACS 64 . 75 . Gh - Phase separation and segregation in model systems (hard spheres, Lennard- 
Jones, etc.) 

Abstract -Molecular dynamics simulations are used to study the phase behavior of a single linear 
multiblock copolymer with blocks of A- and B-type monomers under poor solvent conditions, 
varying the block length N, number of blocks n, and the solvent quality (by variation of the 
temperature T). The fraction / of A-type monomers is kept constant and equal to 0.5, and 
always the lengths of A and B blocks were equal (Na = Nb = N) , as well as the number of blocks 
(tia — riB = n). We identify the three following regimes where: (i) full microphase separation 
between blocks of different type occurs (all blocks of A-type monomers form a single cluster, while 
all blocks of B-type monomers form another), (ii) full microphase separation is observed with a 
certain probability, and (iii) full microphase separation can not take place. For very high number 
of blocks n and very high N (not accessible to our simulations) further investigation is needed. 
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The phase behavior of block copolymers has recently 
attracted much interest in experimental and theoretical 
studies 146. For melts of multiblock copolymer chains 
the phase diagram almost resembles the phase diagram of 
that with diblock copolymer chains; The lamellar struc- 
ture is expected for the most symmetrical case, an ap- 
proach which is particularly valid in the strong segre- 
gation limit .33|- In the case of infinitely dilute so- 
lutions, it is sufficient to consider the behavior of iso- 
lated chains, where interactions (often of short range) 
within the chain are relevant. The case of a single lin- 
ear multiblock copolymer chain |47H51) is very interest- 
ing even when only two different types (A, B) of blocks 
of the same length are considered (see Ref. [51] for dis- 
cussion). Interestingly, multiblock copolymer chains have 
also close relation to the various toy-models (i.e., the HP 
model |52j ) which try to mimic the behavior of various 
biomacromolecules, which are formed by periodically re- 
peated structural units ("monomers") along their chain, 
in order to understand complicated biological phenomena, 
i.e., protein folding [53], helical structures [5J, etc. 

Multiblock copolymers composed of two different types 
of blocks (A and B) under good solvent conditions form 



coil structures and the chain conformations are essentially 
dictated by the repulsive interactions between the differ- 
ent monomers [47j[48]. In fact, the chemical difference of 
monomers should be kept responsible for an expansion in 
the chain dimensions with respect to the equivalent ho- 
mopolymcr chains (same total number of monomers) un- 
der the same thermodynamic conditions [4"5M50j . Also, 
it is expected that the spherical symmetry of the macro- 
molecule should be distorted adopting a slightly ellipsoidal 
overall formation [4"8H5T| . Nevertheless, the most interest- 
ing behavior is expected when the solvent quality is var- 
ied. Under poor solvent conditions the chain collapses, 
and different scenarios of microphase separation between 
the monomers A and B belonging to different blocks are 
conceivable depending on the block length and the number 
of blocks of the chain at a given temperature [J7] . Iden- 
tification of these different states could be based on the 
analysis of the formed clusters as monomers of the same, 
or different, blocks join together. For rather short chains 
(accessible to simulations), simulation techniques provide 
the best way to the understanding of the phase behav- 
ior of a single multiblock copolymer, while a theoretical 
treatment could be an elaborate task. To the best of our 
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Fig. 1: (Color online) Definition of structural parameters de- 
scribing our linear multiblock copolymer chains, n (in this case 
n = 4) is the number of different blocks A and B denoted with 
different color (or grey tone) and TV is the length of each block. 
All the blocks, irrespective of whether they are of type A or B, 
have the same length N. Then the total length of the chains 
is nN. 

knowledge, the phase behavior of a single linear multiblock 
copolymer has not yet been discussed, even for the most 
symmetrical case. 

In the following, we fill this gap and describe large-scale 
molecular dynamics (MD) simulations of an off-lattice 
bead-spring model of a single linear multiblock copoly- 
mer under poor solvent conditions to provide a picture for 
the phase behavior of such macromolccules. We remind 
the reader that one expects that for such a simulated sys- 
tem sharp phase boundaries between different states do 
not exist; rather smooth crossovers are expected. 

In this study, we discuss results for the phase behavior 
of a single symmetric multiblock copolymer with regular 
succession of two different types of blocks (A and B) com- 
posing the linear macromolecule. By using the term "sym- 
metric" we mean that the length of all blocks (A or B) was 
always equal (Na = Nb = N) and the total number of 
A (ua) and B (ng) blocks was n (n = tia + ns), with 
riA = riB in all cases. Therefore, the fraction of monomers 
of type A and B was always constant and equal to / = 0.5. 
A schematic picture of our parameters is given in fig. [T] 

Then, our chains are modelled by the standard bead- 
spring model }55H58| , where all beads interact with a trun- 
cated and shifted Lennard- Jones (LJ) potential 

Uf J (r)=4ef J [(af J /rr-(af J /r) 6 } + C, r<r c , 

(1) 

where a, (3 = A, B denote the different type of monomers, 
and the constant C is defined such that the potential is 
continuous at the cut-off (r c = 2.5) [60]. For simplicity, 
afj = 1, k B = 1, but efj = eff = 2eff = 1, in order 
to create an unmixing tendency between monomers A and 
B belonging to different blocks as done in previous stud- 
ies Therefore, Ae = eff - l/2(ef^ + eff) was 
kept the same throughout our simulations. The connectiv- 
ity along the chain is maintained by the "finite extensible 
non-linear elastic" (FENE) potential 

t^FENE = --fcrg ln[l - (r/r ) 2 ], < r < r . (2) 

Ufene(t > t"q) = oo, and the standard choices ro = 1.5, 
k = 30 were used [55H62]. 

For this model, the O temperature is known only rather 
roughly, namely O = 3.0 Being interested in T < O, 
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Fig. 2: (Color online) Snapshot pictures of three different 
multiblock copolymers of the same total length nN — 600 at a 
low temperature T = 1.5 ((a): N = 6, (b): TV = 15, and (c): 
N = 60). In case (c) we have the formation for two clusters 
of different blocks which are always phase separated. In case 
(a), full phase separation does not take place and the number of 
clusters N c i has a symmetric variation around an average value 
2 < N c i < n. The case (b) is the intermediate case, where full 
microphase separation occurs with a probability P(N c i) > 
(N c i = 2 / n). Different colors (or grey tone) correspond to 
different type of monomers (A,B). 

we simulate temperatures ranging from 1.5 to 3.0, where 
the chain collapses, and monomers of A and B blocks clus- 
ter together with monomers of the same type and mi- 
crophase separation between A and B monomers takes 
place. Using standard MD simulations [5"5H6"T| . the tem- 
perature was controlled by the Langevin thermostat, i.e., 
the equation of motion 

d 2 fi dfi - . 

mLj-Tp = -vUi-rriLji— (3) 

for the time dependence of the monomer coordinates r, (t) 
is numerically solved. Here t is the time, Ui is the to- 
tal potential the i-th bead experiences and the friction 
coefficient 7 was chosen 7 = 0.5 (tolj is also taken as 
unity for simplicity). The white noise Ti(t) satisfies the 
standard fluctuation dissipation relation, (Ti(t) ■ Tj(t')) = 
6kBTjSijS(t — t'). To equilibrate the systems close to the 
temperature, simulations of typically 30 x 10 6 r were car- 
ried out at T = 3.0, using a leap frog algorithm |63j with 
an integration time step of At = 0.006r, where r is the 
natural time unit. After equilibration, we collect a num- 
ber of independent samples (typically 500), which we use 
as initial configurations for slow cooling runs. For longer 
chains, temperatures higher than T = 3.0 were used in 
order to facilitate the procedure of obtaining initial inde- 
pendent samples. We note here that the solvent is taken 
into account in our model only implicitly by tuning the 
temperature as done in previous work (5"5Tl61j . Then, de- 
crease of the temperature corresponds to higher incom- 
patibility of the implicit solvent with the monomers. For 
each cooling history, we lower the temperature in steps of 
AT = 0.1 and simulate the system at each temperature 
for a period that exceeds the relaxation time of the chains, 
using the final configuration at each (higher) temperature 
as starting configuration for the next (lower) temperature. 
At low enough temperatures (typically T < 2.2), where 
dense "clusters" of a few neighboring blocks are formed, it 
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Fig. 3: (Color online) The probability of formation of a single 
cluster of monomers A and a single cluster of monomers B in a 
linear symmetric multiblock copolymer (P(N c i), N c i — 2 7^ n) 
is plotted for different values of block length N, number of 
blocks n and temperature T. The different shading of colors 
corresponds to different values of this probability. 

is not possible to run simple MD simulations long enough 
to sample the phase space adequately. Therefore, using 
this procedure of independent cooling histories is indis- 
pensable for obtaining meaningful statistical results. This 
large statistical effort prevents the study of exceedingly 
long chains, where also the relaxation time of the chains 
becomes very high. We have also performed runs where 
the temperature was changed in one step from T = 3.0 to 
the desired (lower) temperature, and after equilibration, 
properties were calculated. While for T > 2.0 the resulting 
properties did not show any significant dependence on the 
different "history of sample preparation" , somewhat more 
pronounced correlations are seen for temperatures lower 
than T = 2.0 for a single run, which is an indication that 
the intrinsic relaxation times begin to exceed the available 
simulation time scale. Therefore, the use of independent 
cooling histories was necessary in order to perform a mean- 
ingful statistical analysis down to temperatures T = 1.5. 

When the chain collapses at the low temperatures 
and microphase separated clusters of different blocks are 
formed, in order to define the different types of phase be- 
havior one needs to identify for each configuration that is 
analyzed, which blocks belong to a cluster. As in previous 
studies [64] , we have used the standard Stillinger neighbor- 
hood criterion for monomers |65| . When two monomers 
are less than a distance r« apart, they belong to the same 
cluster, but monomers A and B are always considered to 
belong to different clusters. We followed the standard 
choice r n = 1.5, but we have also checked smaller values 
for r n (i.e., r n = 1.2) and similar qualitative results were 
obtained. Higher values of r n are physically hardly signif- 
icant, since then the monomers are too weakly bound, due 
to the rapid fall-off of the LJ potential, fig. [5] shows the 
different types of microphase separation that can be ob- 



Fig. 4: (Color online) Same as in fig. [3] but the colors are 
considered only according to the four cases for the probability 
P(Nci), (N c i — 2 7^ n) as indicated for clarity. 



served in the case of symmetric multiblock copolymers. In 
case (a) , blocks of the same type can join together form- 
ing clusters, but never all the blocks of the same type 
are able to form a single cluster phase separated from an- 
other cluster of B monomers. In this case, the number 
of formed clusters N c i fluctuates symmetrically around an 
average value (2 < N c i < n). This behavior is also an in- 
dication that the adopted simulation protocol that we fol- 
lowed in this work leads to a meaningful statistical analysis 
by considering this large number of independent cooling 
histories, as is done in previous work 56,59,61.64 . More- 
over, analyzing separately the properties for clusters of A 
monomers and for clusters of B monomers we obtain the 
same results, which also shows the validity of our simula- 
tion procedure. In case (c), all monomers A of all different 
blocks belong to a single cluster of A-type monomers and 
another cluster of B-type monomers is also formed. Then 
case (b) is the intermediate case, where strong fluctua- 
tions in the number of clusters is observed and the occur- 
rence of full phase separation takes place with a proba- 
bility P(N cl ) > 0.0, N d = 2. Thus, the different phase 
behavior in multiblock copolymers can be characterized 
by defining the boundaries between cases (a), (b), and (c) 
in terms of the probability P(N c i) for N c i = 2. 

figs. [3] and H] present results of the probability of de- 
tecting full phase separation (i.e., all blocks of the A-type 
monomers form a single cluster, while all monomers of B- 
type belong to another cluster (P(N c i), N c i = 2) between 
monomers A and monomers B. The shading of the colors 
in fig. |3] corresponds to this probability, whereas in fig. H] 
only the colors of the four different cases are shown for 
clarity, without color shading according to the probability 
P(N c i = 2), as indicated in the figures. Considering this 
probability for the full phase separation, we are able to 
describe the different cases of the phase behavior, which 
fig. [2] illustrates with characteristic snapshot pictures. We 
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have considered a broad variation of parameters, n, N, 
and T, with T < 0, and we show here representative re- 
sults of our simulations for the description of the different 
regimes. We can observe that for small block length N 
(N < 20), full phase separation (P(N c i = 2) = 1, as in 
case (a) of fig. can be hardly seen independently of 
the number of blocks n, even at low temperatures, where 
the chains collapse and form globules. In this regime, the 
probability of full phase separation is P(N c i = 2) < 0.5. 
Looking more carefully in fig. [3] we can actually observe 
that decrease of n increases the probability of the forma- 
tion of single clusters (one with A monomers and another 
with B monomers), as expected. 

The most interesting regime, where full microphase sep- 
aration occurs, is for rather small values of n, 4 < n < 20 
(of course for n — 2 the clusters of A and B monomers are 
considered always separated) . For this to happen, also the 
length of the blocks N has to be higher than a certain value 
(N > 20). This regime is indicated in figs. [3] and fig. S] 
with yellow color (or the very light grey tone). It turns 
out that the increase of N for certain value of n favors full 
phase separation of A and B blocks, since for higher N 
full separation is observed, even at higher temperatures. 
In particular, for n = 4,8, 12 and N = 150 such an effect 
takes place even at temperatures in the range T = 2.3 to 
T = 2.4. At high temperatures, independently of n and 
N, the probability P(N cl ) « (N cl = 2). In this case, the 
different blocks prefer to be apart at high temperatures 
and can form clusters with other blocks with monomers of 
the same type only occasionally, due to thermal fluctua- 
tions. To conclude, full phase separation is favored when 
the number of blocks is low, and the block length higher 
than a certain value (N sa 20), for the range of tempera- 
tures that we used in our simulations. We note that the 
different states occur gradually and the boundaries be- 
tween the different states are to be considered as smooth 
crossovers (indicated also by the results of fig. [3]). How- 
ever, we have observed with closer inspection of our data 
that these crossovers become sharper with the increase of 
the block length N. 

Very interesting is also the phase behavior observed for 
high number of blocks n, while N remains finite. In this 
case, n is so high that the occurrence of full phase sepa- 
ration, as we have defined it already in our discussion, is 
rather not observed for the given set of potential parame- 
ters (Ae). What actually one sees is domains of different 
(A and B) microphase separated clusters for the range 
of N presented in this study. It would be very interest- 
ing to simulate very long block lengths, in order to check 
if full microphase separation is possible for very high iV 
(or changing x ~ Ae/T; changing cither x or N would 
result in tremendous difficulties in simulating such a sys- 
tem) keeping also high the number of blocks n. How- 
ever, if one naively tries to extrapolate our data for n and 
N in the range where full microphase separation occurs, 
it could be suggested that full phase separation for high 
n could require also very high values of N. It has also 



been discussed for copolymer melts [1] , that the geometry 
of the microphase separated regions is controlled by the 
number of blocks n, as well as other parameters, i.e., rela- 
tive size and arrangement of the blocks. Nevertheless, one 
would expect in the long chain limit in equilibrium that 
a ground-state type structure would be a single lamellar 
domain, where an interface between all A-type and B-type 
blocks is formed, similarly to what is known for multiblock 
copolymer melts [33]. Such a structure would have much 
less (unfavorable) A-B contacts rather than a multidomain 
structure of A and B clusters, which is kinetically favored 
due to the preparation protocol in the simulation and the 
chosen set of parameters accessible to simulations. 

In summary, we have discussed the microphase sepa- 
ration of symmetric linear multiblock copolymers, where 
a chain contains blocks of monomers of types A and B, 
which alternate along the chain via MD simulations of an 
off-lattice bead spring model. The fraction of monomers 
A (Ja) and monomers B (fs) was always constant at 
/ = /a = /b = 0.5. Varying the number of blocks n, the 
block length N, and the temperature T, we were able to 
detect the regimes where full phase separation (all blocks 
of type A belong to a single cluster and all blocks of type B 
to another cluster), as well the regime that full phase sep- 
aration does not take place. We found that rather small 
values of n (4 < n < 20) and increasing block length 
N (N > 20) favors full separation of the blocks, but N 
should in no case be lower than the aforementioned value. 
Let us point out here that, the behavior of multiblock 
copolymers can be parallelized with that of various bio- 
logical macromolecules, which are formed by periodically 
repeated structural units ("monomers") along the chain. 
In such systems the capability of the formation of clus- 
ters could be related to the probability of forming bonds 
and occurrence of interactions between different blocks on 
the same chain. Further investigation is needed for very 
high number of blocks and very high block lengths, not 
accessible to our simulations. 
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